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ABSTRACT 

We present a method for computing the probability distribution of microlensed 
light curve derivatives both in the case of a static lens with a transverse velocity, and 
in the case of microlensing that is produced through stellar proper motions. The distri- 
butions are closely related in form, and can be considered equivalent after appropriate 
scaling of the input transverse velocity. The comparison of the distributions in this 
manner provides a consistent way to consider the relative contribution to microlens- 
ing (both large and small fluctuations) of the two classes of motion, a problem that 
is otherwise an extremely expensive computational exercise. We find that the relative 
contribution of stellar proper motions to the microlensing rate is independent of the 
mass function assumed for the microlenses, but is a function of optical depth and shear. 
We find that stellar proper motions produce a higher overall microlensing rate than a 
transverse velocity of the same magnitude. This effect becomes more pronounced at 
higher optical depth. With the introduction of shear, the relative rates of microlensing 
become dependent on the direction of the transverse velocity. This may have impor- 
tant consequences in the case of quadruply lensed quasars such as Q2237+0305, where 
the alignment of the shear vector with the source trajectory varies between images. 

Key words: gravitational lensing - microlensing - numerical methods. 



1 INTRODUCTION 

Calculations of the properties of a light curve from a grav- 
itationally microlensed quasar have traditionally assumed 
that the variation in observed flux is the result of an effec- 
tive transverse velocity of the source with respect to the 
foreground galactic starfield. These calculations have as- 
sumed that the contribution to microlensing from the ran- 
dom proper motions of stars is a negligible one because the 
galactic transverse velocity is an order of magnitude larger 
than the typical stellar proper motion. However in the case 
of Q2237+0305, the only object in which cosmological mi- 
crolensing has been confirmed (Irwin et al. 1989; Corrigan 
et al. 1991), the line-of-sight velocity dispersion is compa- 
rable to the likely transverse velocity. Foltz et al. (1992) 
measured the central velocity dispersion of Q2237+0305 to 
be ~ 215fcmsec _1 . This is to be compared with the esti- 
mated galactic transverse velocity of ~ 600 km sec -1 , and 
the suggestion by Witt & Mao (1993) who provide evidence 
from the observed high magnification events (HMEs) for a 



smaller transverse velocity. Q2237+0305 is the perfect ob- 
ject from which to study microlensing. The source quasar is 
lensed by an intervening galaxy, producing 4 images, each of 
which is observed through the galactic bulge with an optical 
depth in stars that is of order unity (eg. Schmidt, Webster 
& Lewis 1998). In addition, the close proximity of the lens- 
ing galaxy means that the effective transverse velocity may 
be high. However the potential contribution of the velocity 
dispersion makes accurate direct modelling of the microlens- 
ing light curves difficult. Therefore while Q2237+0305 is the 
best object from which to observe microlensing, a realistic 
numerical analysis of the observations is difficult. 

This paper describes the relative contributions to mi- 
crolensing of random stellar proper motions, and of a trans- 
verse velocity of the entire lens. Two approaches have been 
used to tackle this problem in the past. Kundic & Wambs- 
ganss (1993), and Wambsganss & Kundic (1995) computed 
light curves that included the effect of stellar proper motions 
by producing a magnification map at many successive inter- 
vals for an evolving field of point masses. Light curves that 



© 1998 RAS 



2 J. b. t>. Wyithe, K. L. Webster & h. L. lurner 



resulted from the random motion of stars alone were then 
produced by looking at the changing magnification at a fixed 
source point, and compared with light curves that resulted 
from a transverse velocity of the source with respect to one of 
the starfields. These light curves were produced in the usual 
way by shifting the magnification pattern about the source 
point. Schramm et al. (1992) obtained an expression for the 
caustic velocities that result from a stream motion of the 
star field. This idea was extended by Kundic, Witt & Chang 
(1993) who calculated the caustic velocity resulting from a 
stellar velocity dispersion. They calculated the area swept 
out per unit time by the caustic network in both the cases of 
a lens composed of a static field of point masses with a trans- 
verse velocity, and one where the the point masses are given 
a random proper motion. These papers describe the effect 
that proper motions have on the statistics of high magnifi- 
cation events (HMEs). Unfortunately their conclusions are 
not consistent. Kundic, Witt & Chang (1993) claim that 
the relative efficiency of stellar proper motion and trans- 
verse velocity in terms of the resulting frequency of HMEs 
is a function of the optical depth, while Kundic & Wambs- 
ganss (1993) find no such dependence. The effect of the in- 
clusion of a shear in the models is mentioned in these papers. 
However no conclusions are presented regarding variation in 
the relative importance of microlensing due to stellar proper 
motions with the direction of the transverse velocity. 

Our approach is to compute the distribution of light 
curve derivatives in the cases of a static lens with a trans- 
verse velocity, and of a static source point in a model where 
the point masses are given a random proper motion. As 
shown below these distributions are very similar up to a 
scaling factor in the derivative and can be quantitatively 
compared in order to find the relative contributions to mi- 
crolensing of the two cases. The advantages of this approach 
are two fold. Firstly, the average distribution of light-curve 
derivatives can be computed from points that are not se- 
quential along a light curve. The distribution of light-curve 
derivatives for the cases of proper motion and transverse ve- 
locity can therefore be computed at the same computational 
expense. This is in contrast to the direct computation of 
a light curve (Kundic & Wambsganss 1993; Wambsganss & 
Kundic 1995) where the increase in time is equal to the num- 
ber of points on the light curve. Secondly, the distribution 
of light-curve derivatives provides information on the mi- 
crolensing rate in all derivative regimes. It is therefore more 
suited to the analysis of monitoring data than a method 
that considers only the frequency of caustic crossings be- 
cause HMEs are rare, and may be missed by the low sam- 
pling rate. 

This paper is concerned with the construction of the 
histogram, the nature of microlensing due to stellar proper 
motions, as well as its dependence on optical depth, shear 
and direction of the source trajectory. It is presented in 8 
parts. Sections |^ and ^| discuss the utility of the distribution 
of microlensing light-curve derivatives, and describe the mi- 
crolensing models used in the paper. Sections ^ and ^ de- 
scribe the construction of the derivative histograms for the 
proper motion and transverse velocity cases, and the method 
for their comparison. In Sections |andg the variation with 
optical depth and shear of the relative contribution of proper 
motions and transverse velocity to the microlensing rate are 
discussed. Section H presents the relative rates for the pa- 



rameters of optical depth and shear corresponding to the 
images of Q2237+0305. 



2 APPLICATION OF THE DISTRIBUTION OF 
MICROLENSED LIGHT-CURVE 
DERIVATIVES 

This paper is part of an analysis of the 12 years of pub- 
lished monitoring data for Q2237+0305. In this analysis the 
distribution of light curve derivatives is considered, an ap- 
proach that yields information on both the effective galac- 
tic transverse velocity, the source size and the microlens 
mass function. These aspects are discussed in subsequent pa- 
pers (Wyithe, Webster & Turner 1999b, c; Wyithe, Webster, 
Turner & Mortlock 1999), but the distribution of the deriva- 
tives due to proper motion of stars is an integral part of the 
analysis. The gain in computational speed makes consider- 
ation of proper motions in a large range of models feasible. 
However, the histogram of microlensed light curve deriva- 
tives does not contain information on the structure of the 
light-curve. Therefore, to apply the derivative distributions 
due to proper motions to a set of monitoring data we make 
the assumption that a typical section of light-curve result- 
ing from stellar proper motions will have similar statistics 
to a section of light curve produced by a static lens in com- 
bination with an appropriate transverse velocity. Further- 
more, for the distribution of derivatives to be applicable to 
short sections of light-curve, caustic crossings in the two 
cases should have the same clustering properties. 

To justify these assertions we note firstly that (as shown 
below) at an appropriate transverse velocity, light curves 
from a static field have a very similar average distribution 
of microlensed light-curve derivatives to that produced by 
proper motions, indicating equal average rates of light curve 
peaks and average peak durations (the peak heights must 
also be equal). In addition, the similarity of the two classes 
of derivative distribution suggests an approximately equal 
temporal clustering of peaks since additional clustering in 
one case would produce an excess of small derivatives and a 
corresponding dearth of mid-range derivatives over the other 
case. Secondly, at any given time caustics that move under 
the influence of random proper motions have the same spa- 
tial clustering since there is no correlation between motions 
of individual point masses. Also, in general the motion of 
two caustics that are in close proximity are governed pri- 
marily by the motions of separate stars or groups of stars. 
Caustics that have spatial correlation therefore do not have 
a corresponding correlation in the direction or magnitude 
of their velocity, and hence there should be no additional 
temporal correlation of caustic crossings. 

With these points in mind, we subsequently assume that 
the temporal clustering of caustic crossings is the same for 
the two classes of motion, and that light-curves produced 
by proper motion and by an appropriate transverse velocity 
can be considered equivalent. Sampled model light-curves 
due to stellar proper motions can then be produced from 
static model light curves with the appropriate transverse 
velocity and compared with observation. 
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Table 1. The parameters of microlensing models considered in this paper. 



simulation no. 


K 


7 = (71,72) 


simulation type 


source line length 


no. stars 


no. simulations 




(Hth) 


1 


0.025 


(0,0) 


MSOLAR 


50r; o 


500 


500 


1.04±.01 


1.05 


2 


0.05 


(0,0) 


MSOLAR 


50>7 o 


500 


500 


1.08±.01 


1.11 


3 


0.1 


(0,0) 


MSOLAR 


50r; o 


500 


500 


1.18±.01 


1.23 


4 


0.2 


(0,0) 


MSOLAR 


50»7 


511 


250 


1.54±.04 


1.56 


5 


0.3 


(0,0) 


MSOLAR 


50r; o 


1127 


250 


2.02±.05 


2.04 


6 


0.1 


(0,0) 


MSOLAR 


25r/ D 


1071 


250 


2.68±.05 


2.78 


7 


0.5 


(0,0) 


MSOLAR 


25r/ 


2194 


250 


4.09±.10 


4.00 


8 


0.6 


(0,0) 


MSOLAR 


25r/ D 


4597 


100 


6.16±.15 


6.25 


9 


0.2 


(0,0) 


SALPETER 


50r/ o 


1792 


200 


1.51±.02 


1.56 


10 


0.3 


(0,0) 


SALPETER 


50r/ o 


3792 


200 


2.03±.05 


2.04 


11 


0,1 


(0,0) 


SALPETER 


50?? o 


7333 


200 


2.63±.05 


2.78 


12 


0.5 


(0,0) 


SALPETER 


25i7o 


5611 


100 


3.95±.12 


4.00 


13 


0.6 


(0,0) 


SALPETER 


10r/ o 


5423 


100 


6.26 ±.33 


6.25 


11 


0.3 


(-0.4,0) 


MSOLAR 


50r/ o 


855 


200 


2.95±.09 


3.03 


15 


0.3 


(-0.2,0) 


MSOLAR 


50r/ o 


789 


200 


2.18±.04 


2.22 


16 


0.3 


(0.2,0) 


MSOLAR 


50?? o 


2102 


200 


2.23±.03 


2.22 


IV 


0.3 


(0.4,0) 


MSOLAR 


50»? o 


5739 


200 


3.00±.01 


3.03 


18 


0.36 


(0.4,0) 


SALPETER 


10»? o 


2926 


60 


3.92±.08 


4.01 


19 


0.36 


(0.28,0.28) 


SALPETER 


10r/ o 


1774 


60 


3.95±.14 


4.01 


20 


0.36 


(-0.4,0) 


SALPETER 


10r, o 


643 


60 


3.75±.15 


4.01 


21 


0.69 


(0.71,0) 


SALPETER 


10r/ o 


2782 


30 


2.52±.17 


2.45 


22 


0.69 


(0.50,0.50) 


SALPETER 


10r7 o 


2311 


30 


2.46±.10 


2.45 


23 


0.69 


(-0.71,0) 


SALPETER 


10»? o 


1058 


30 


2.42±.10 


2.45 


24 


0.59 


(0.61,0) 


SALPETER 


10r/ o 


8583 


30 


4.77±.10 


4.90 


25 


0.59 


(0.43,0.43) 


SALPETER 


10r, o 


6675 


30 


4.82±.29 


4.90 


2(5 


0.59 


(-0.61,0) 


SALPETER 




2219 


30 


4.62±.44 


4.90 



3 THE MICROLENSING MODELS 

Throughout the paper, standard notation for gravitational 
lensing is used. The Einstein radius of a 1M@ star in the 
source and image planes are denoted by rj and £ respec- 
tively. The normalised shear due to external mass is denoted 
by 7 — (71,72), and the convergence or optical depth by n. 
The model for gravitational microlensing consists of a very 
large sheet of point masses that simulates the section of 
galaxy along the image line of sight, together with a shear 
term that includes the perturbing effect of the mass distri- 
bution of the lensing galaxy as a whole. The normalised lens 
equation for a field of point masses with an applied shear in 
terms of these quantities is 



V = 



- 71 
-72 



-72 
f 71 



x + 



3=0 



(x j 



(1) 



Here x and y are the normalised image and source positions 
respectively, and the xl are the normalised positions of the 
point masses. The magnification, or increase in flux due to 
gravitational lensing of a point source is obtained from the 
inverse of the Jacobian determinant of Eqn |l|: 

1 



Hp 



det(A(x)) ' 



(2) 



where A = is the Jacobian matrix. The magnification 

Hp may be negative, corresponding to the magnification of 
an inverted image. In the case of multiple imaging of a point 
source the magnification Hv is given by 



Hp 



^ \f-Limages 



(3) 



We consider two classes of model in our analysis, both 
consisting of point masses with positions that are distributed 
randomly in a disc. In sections ^| and |^ each model star is as- 
signed a mass of IMq . Thes e simula tions are termed to be of 
type MSOLAR. In sections 6^ , sl}|andp|the model is based 
on that of Wambsganss, Paczynski & Katz (1989), in which 
masses are distributed according to a Salpeter mass function 
(p(m)dm oc m -2 ' 35 ) in the range 0.1M Q < mi < 1.0M Q . 
The label SALPETER refers to simulations of this type. 
None of the models presented in this paper include a compo- 
nent of continuously distributed matter k c in addition to an 
optical depth in stars k, . The results can however be applied 
to this case since a model that includes continuous matter is 
mathematically equivalent (through the parameter transfor- 
mation of Paczynski (1986)) to a model that contains only 
compact objects, but which is described by a different set 
of microlensing parameters («' and 7') and different scale 
lengths (£0 and jj' ). For the case where k c < 1: 



>lo 



1 — K c 

7 x (1 

£ 



rj X s/l — k c 



The point source magnifications are related by 
// = (! — k c ) 2 h- 



(4) 



(5) 



The region of the lens plane in which image solutions 
need to be found to ensure that 99% of the total macro- 
image flux is recovered from a source point was described 
by Katz, Balbus & Paczynski (1986). In the presence of an 
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applied shear these regions are elliptical in shape. The union 
of the areas of the lens plane that correspond to the flux col- 
lection area of each point on the source line is known as the 
shooting region. The method of determining the dimensions 
of the shooting region is described in Lewis & Irwin (1995), 
and Wyithe & Webster (1998). The radius of the disc of 
point masses is chosen to be 1.2x that required to cover 
this shooting region. Eqn [l] is solved through the inversion 
technique of Lewis et al. (1993) and Witt (1993). This tech- 
nique finds all image solutions of a source line that lie within 
the shooting region using the fact that these lie on one of 
two classes of image curve. The first of these curves is a line 
that tends asymptotically to the unperturbed (by the point 
masses) image of the source line far from the starfields edge. 
This line may pass over one or more of the stars. The second 
class of image curve is a collection of loops, each of which 
passes over one or more stars. An example of the lensed 
images of a line is shown in Figure jjj. We implement the 
contouring method as described by Lewis et al. (1993). Im- 
age solutions are found iteratively along the image curves. 
When each new image point is located the image magnifi- 
cations are added to predefined source points through an 
interpolation between the current and previous steps. The 
interpolation schemes between two source points (i and 
were given by Lewis et al. (1993): 



y OC Hivtiage.1 



where 



= — 2 for fit or > 



(1-k) 2 - 7 2 



Q 

a 



for Hi and < 0.2, 



= 1 otherwise. 



(6) 



(7) 

(8) 
(9) 



Table ^ displays the parameters of simulations discussed 
in this paper. Each simulation is labelled numerically, and 
is described by the associated values of k and 7, as well as 
the simulation type. Columns 5 and 6 in Table |l| show the 
length of source track and the number of stars in each indi- 
vidual simulation respectively. Column 7 shows the number 
of individual simulations in the set. The final 2 columns 
show the mean amplification resulting from the set of simu- 
lations, and the theoretical value (jlth = |(1 — k) 2 — 7| 2 | _1 ) 
for comparison. Each set of simulations was divided into 5 
subsets, and the error in the mean magnification of the mod- 
els estimated from the standard deviation in the resulting 
values. Figure ^ demonstrates that the magnification distri- 
butions of our simulations are independent of mass function 
and source trajectory direction as required. 



4 CONSTRUCTION OF THE HISTOGRAM OF 
DERIVATIVES 

To calculate the distribution of light curve derivatives that 
result from stellar proper motions, the derivative of the mag- 
nification is taken with respect to time at the position of 
each image. These derivatives are then added to give the 
total change in magnification due to the movement of the 
point masses. Eqn hi can be written in component form: 



2/i 



and 



(1-7)2:1 - 



Xi — xi 



3=0 



(xi — x{) 2 + (X2 — X ] 2 ) 2 



(1 + j)x 2 - m? 



3=0 



X2 — X 2 
(Xl — x{) 2 + (X2 — x{) 2 ' 



(10) 



(11) 



where x\ , m? ,Xi,yi, for i = 1, 2 are the position and mass of 
the jth star, and the position of the image and source respec- 
tively. For simplicity the source and image axes have been 
aligned with the shear vector in this expression (7 = (7, 0)) 
in order to diagonalise the shear matrix. The positions of 
the point masses are functions of time: 

dr j 

vm^. (12) 

For a given image we write down the expression for the 
derivative of magnification with respect to time in terms 
of these velocities: 



~dt 



d Hp 



+ 



+ 



dx± dx{ sr--\ dxi dx J 2 

~dx~i \^JxJ~^ + ^'dxJ~dT 

\i=o 1 j=o 2 

d Hp ( d X2 dx{ 8x2 dx\ 

~dx~2 \^~d^ 1 ~dT + l^J^~dt 

\j=o 1 3=0 2 , 



Ed H P dx\ \—\ d Hp dx 2 
~dxJ~dT + ^ Qxi, dt 

3=0 1 3=0 2 



(13) 



To compute the derivatives for i, k = 1,2 the derivatives 

8 x k 

Pi- = and 4^2- = are taken simultaneously. The full 



expression for -jf-) pr oper is given in appendix (Eqn |A§) 



As with the magnifications, the light curve derivatives 
must be interpolated between iterations onto predefined 
source points. Where the magnification derivative is the re- 
sult of stellar proper motions, the interpolation is found by 
noting that if y tx Hp then y oc (h p + ^f-At) a also holds for 
small At. Following this we obtain an interpolation scheme 
for the derivative: 



V oc Hp 



d Hp 
~~dT' 



(14) 



From Eqn ^, the time derivative of the point source magni- 
fication is 



d-Hv • / ^ d Hi 

— = ^ si gn(H i )^ T . 

i=0 



(15) 



The magnification is not directly measurable, however we 
can find the magnification derivative in terms of the resul- 
tant rate of change of magnitude. 



(16) 



d Am 1.09 d/ip 
dt Hp dt 

The histogram of derivatives for a static microlensing 
model can be calculated in a straight forward manner di- 
rectly from the combination of a computed light curve with 
a transverse velocity of the lens (v tra n = ^f). However, in 
order to obtain a de rivat ive of comparable accuracy to that 
calculated from Eqn A2, we find the rate of change in mag- 



nification in this case analytically as before. 
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Source position (Einstein radii) 

Figure 1. The images of a source line (top) and point source 
light curve (bottom) for a simulation having an optical depth of 
k = 0.4. The source trajectory is aligned at 45 degrees to the 
shear vector. The dashed light line shows the theoretical average 
magnification for this model. 



Figure 2. The cumulative histogram of derivatives (in units of 
magnitudes per year) calculated from the light curve in Figure 
|l] (dot-dashed line). Each of the point masses in this simulation 
were given a velocity parallel, opposite in direction, and of the 
same magnitude as the transverse velocity. The distribution of 
derivatives wa s th en calculated for this same simulation using 
the expression A2, and is shown as the solid light line 



dt J tran 



dfipdx! d [ip dx 2 \ dyi 
dx± dyi dx2 dyi J dt 



(17) 



In this expression , and |^ are evaluated by taking the 
derivatives 4^ = 1 and 4-^- = m simultaneously where the 

dyi dyx ■> 

source track is described by y 2 = myi. The full expression 
for -jf-)tran is given in appendix f\] (Eqn |A|). The interpo- 
lation scheme in this case is simply the time derivative of 
that for the computed light curve and so Eqn [l^ also holds 
in this case. In addition, Eqns [TB] and |l^ are correct in this 
case. 

Kundic, Witt &i Chang (1993) describe the relation- 
ship between the caustic velocity that results from a galactic 
transverse velocity (vtran), and that resulting from a stream 
velocity of the local starfield {v stream) - The microlensed light 
curve derivative is proportional to the caustic velocity, and 
so we find for a static starfield having both a transverse and 
a stream motion that 



dt 



Vfr 



+ 



1-7 






1+7 



lustre 



(18) 



In this work, we restrict our attention to the investigation of 
the relationship between the microlensing rate due to proper 
motions and galactic transverse velocity. Note however that 
the results can easily be converted to describe the corre- 
sponding relationship with a stream motion. 



5 THE METHOD OF COMPARISON 

We investigate the contribution of stellar proper motions to 
microlensing statistics by looking for the transverse velocity 
(equivalent transverse velocity) that produces a derivative 
histogram closest to the one produced by the proper motions 
of the point masses in the same fields. The similarity of two 
cumulative distributions Pi and P 2 is quantified by their KS 
difference D where 



D = max(\P± - P 2 \). 



(19) 



This procedure leads to the natural definition of a constant, 
which following Kundic, Witt & Chang (1993), and Wamb- 
sganss & Kundic (1993) we call the effectiveness parameter 

_ equivalent transverse velocity 
tot velocity dispersion 

Note that this value differs to that defined by Wambsganss 
& Kundic (1993), as it takes account of all microlensing as 
opposed to only HMEs. For the sake of clarity we therefore 
refer to the quantity defined by Eqn ^ as the total effec- 
tiveness parameter (a to t)- 

Figure ^ shows the images of a source line as well as 
the corresponding point source light curve for a simulation 
having an optical depth of k = 0.4. The source trajectory is 
aligned at 45 degrees to the zi-axis. Figure [| displays the 
cumulative histogram of derivatives that was calculated di- 
rectly from the combination of this light curve with a trans- 
verse velocity (dot-dashed line) . Each of the point masses in 
this simulation were given a velocity equal to the transverse 
velocity and parallel to the source line. In the absence of 
a shear, the transverse and stream motions are equivalent. 
The distribution of derivatives at stationary source points 
alon g th e source line was then calculated using the expres- 
sion A2, and is shown as the solid light line in Figure H. 
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In addition, the histogram of derivatives in the transverse 
motion case was also calculated using the combination of 
the analytical expression for the light curve derivative (Eqn 
A3), with the transverse velocity. 

The histograms calculated from the analytical deriva- 
tives agree to within the numerical resolution of the com- 
puter (as they must if our expressions are correct). The 
histograms calculated from the analytical and numerical 
derivatives (calculated from a 5-point derivative along the 
light-curve) display excellent agreement, having a maximum 
KS difference of 0.008. Figure H shows the simple nature of 
the derivative distribution, even when computed over a rel- 
atively short sample length. This is in stark contrast to the 
highly disordered parent light curve. 

As a further check of our method, the stars are again 
given parallel motions resulting in a stream velocity, but this 
time the direction of this motion is perpendicular to that of 
the transverse velocity. In this case our simulation consists of 
100 fields of point masses, each with an optical depth of 0.4. 
The best fit was obtained between the histogram of deriva- 
tives due to stream motions, and that due to a transverse 
velocity, by finding the equivalent transverse velocity that 
minimised the KS difference between the histograms. This 
procedure produced a ratio between the value of the stel- 
lar stream motion and the calculated equivalent transverse 
motion of 0.9993 with a minimised KS difference of 0.007. 
These simulations demonstrate that the method employed 
is working correctly. 



6 VARIATION WITH OPTICAL DEPTH 
6.1 The effect of optical depth 

In order to investigate the dependence of the relative ef- 
fects of a transverse motion and that of a collection of ran- 
dom proper motions, we start by following the example of 
Kundic, Witt & Chang (1993) who explore the unrealistic 
situation in which the point masses are all solar mass objects 
that have a motion of fixed magnitude, but a random direc- 
tion lying in the lens plane. We label simulations with this 
form for the stellar velocity distribution by CONST. DISP. 
These models allow us to explore the relationship between 
the microlensing properties of the two types of motion in a 
simple situation without the added complexities of velocity 
and mass functions. These are discussed in sections 3.2 and 



We consider the effect of optical depth on the value 
of a to t- M SOLAR type simulations were made at optical 
depths between k = 0.025 and k — 0.6 (simulations 1-8 in 
Table 0). Table ^ displays the value of atot obtained as well as 
the best fit KS difference, and the derivative at which the KS 
difference was found (in brackets) . Table ^ also displays the 
results of sections 5.2 and 6.3. The quoted error in these, and 



all values of a to t presented in this paper have been estimated 
by dividing the total simulation at each optical depth into 5 
smaller simulations, and calculating the standard deviation 
in the resulting values. 

Figure ^ shows graphically the relationship between k 
and a to t- The results of this section are represented by the 
thin dark line on this plot. We find that at high optical depth 
the effect of random proper motions on the microlensing 




OPTICAL DEPTH 

Figure 3. The relationship between optical depth n and 
the effectiveness parameter atot- In these simulations there 
was no applied shear (7 = 0). The thin, thick and 
dashed lines show the dependence for simulations of types 
MSOLAR/CONST.DISP., SALPETER/CON ST. DISP. and 
SALPETER/GAUSS. DISP. respectively. 



Table 3. Results from simulations of type M SOLAR, with veloc- 
ity dispersions that are of type CONST. DISP. The simulations 
had no applied shear (7 = 0). Details are found in the text. 



Optical 
Depth 



Effectiveness Parameter Effectiveness Parameter 



a HME ( 



d Am 
dt 



> 0.15) 



O-HME ( d jf 1 > 1) 



0.10 
0.20 
0.30 
0.40 
0.50 
0.60 



1.02 
1.09 
1.16 
1.20 
1.30 
1.34 



1.01 
1.08 
1.15 
1.13 
1.22 
1.29 



rate is larger than that of a transverse velocity of the same 
magnitude. However at low optical depth the proper motions 
produce a lower microlensing rate. The total effectiveness 
parameter, or any similar quantity must have a value of 1 
in the limit of zero optical depth, and we find evidence that 
there is a trend towards a total effectiveness parameter of 1 
in this limit. This has not been pursued further, because in 
the case of Q2237+0305 we are interested in optical depths 
of between about k = 0.2 and k = 0.6. 

The reasons for the behaviour of a to t at both low and 
high optical depths are illuminated through the compari- 
son of probability density functions of light curve deriva- 
tives produced both by stellar proper and galactic trans- 
verse motions at various optical depths. Figure ^ displays 
these probability density functions for k = 0.2, 0.4 and 0.6 
in the case where both the magnitudes of the proper motions 
and of the transverse velocity are 300 km sec -1 . In this dia- 
gram the dark lines indicate the probability density function 
for light curve derivatives produced through proper motion, 
while the light lines indicate the function for a static field 
with a transverse velocity. All three cases display an excess 
probability of finding a large light curve derivative in the 
case of microlensing due to the proper motion of stars. The 
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Table 2. Results from simulations of sections 



6.2 
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and [j.j . The simulations had no applied shear (7 = 0). Details are found in the text. 







Effectiveness 






Minimised KS 








Parameter (atot) 






Difference 




Optical 


MSOLAR, 


SALPETER, 


SALPETER, 


MSOLAR, 


SALPETER, 


SALPETER, 


Depth 


CONST. DISP. 


CONST. DISP. 


GAUSS. DISP. 


CONST. DISP. 


CONST. DISP. 


GAUSS.DISP. 


0.025 


0.95± .07 


' 




0.011 (.000015) 






0.050 


0.87± .04 






0.010 (.000078) 






0.10 


0.91± .02 






0.009 (.00067) 






0.20 


n nn_l_ no 

o.yy± .03 


1.01±.03 


l.z0±0.01 


n m n / nn7i A 
0.010 (.00/ 1 ) 


0.00/ (.014) 


0.010 (.013) 


0.30 


1.08± .02 


1.10±.03 


1.28±0.03 


0.012 (.014) 


0.009 (.028) 


0.010 (.037) 


0.40 


1.20± .03 


1.17±.02 


1.46±0.04 


0.018 (.083) 


0.008 (.049) 


0.010 (.050) 


0.50 


1.32± .07 


1.31±.06 


1.60±0.03 


0.009 (.131) 


0.007 (.088) 


0.004 (.076) 


0.60 


1.38± .05 


1.39±.06 


1.77±0.04 


0.007 (.132) 


0.006 (.130) 


0.003 (.127) 




DERIVATIVE (MAGNITUDES/YEAR) 

Figure 4. The probability density functions in both the cases of microlensing due to a transverse velocity (light lines), and due to stellar 
proper motions (dark lines). The probability density functions are shown for the cases of n = 0.2 (solid lines), k = 0.4 (dashed lines), 
and k = 0.6 (dot-dashed lines). The simulations had no applied shear (7 = 0). 



increase in microlensing rate at higher optical depth is also 
illustrated by these functions. 

The excess of large derivatives in light curves that re- 
sult from proper motions is quantified in Table ^. In this case 
the HME effectiveness parameter cihae is defined as the ra- 
tio between the speed of the point masses, and that of the 
transverse velocity required to produce the same probabil- 
ity of finding a light curve derivative above some minimum 
level. The values we have chosen for this minimum level are 
d ^ t m = 0.15, which corresponds to the value used by Wamb- 
sganss & Kundic (1993) scaled by a transverse velocity of 
300 km sec -1 , and d ^ m = 1. This measure differs from that 

' a t 



employed by Wambsganss & Kundic (1993) however because 
we do not require a minimum magnification shift, but rather 
just a minimum rate of change. If these large derivatives are 
only found during HMEs, then an me approximates the ratio 
of probabilities between the two cases that the source is un- 
dergoing an HME at any given time. This is not necessarily 
the same thing as the ratio of the number of HMEs produced 
by the two types of motion because it does not take into ac- 
count the typical HME time-scale. The results presented in 
Table |^ show that the trend of an increasing effectiveness 
parameter with optical depth is present in the high deriva- 
tive regime, consistent with the results of Kundic, Witt & 
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Figure 5. Magnification distributions. Left: k = 0.3, 7 = 0, 
MSOLAR (light line) and SALPETER (dark line). Right: k = 
0.3, 7 = +0.4 (light line) and 7 = -0.4 (dark line). 



Chang (1993). The excess of large derivatives is expected 
following the calculations of Wambsganss & Kundic (1995) 
who show that events produced by random proper motions 
are on average shorter, more frequent and have larger slopes. 

Figure H shows a change in behaviour in the low deriva- 
tive regime. At the higher optical depths considered, the 
increase in large derivatives is offset by a corresponding de- 
crease in the probability of finding a low derivative. However, 
at the lower end of the range there is an increase both in 
the probability of finding the highest, and the lowest deriva- 
tives in the proper motion case. This increased population 
of small derivatives dominates the evaluation of atot at low 
optical depth because the lowest derivatives are by far the 
most common. This explains why the total effectiveness pa- 
rameter dips below 1 at low optical depth. Figure ^| suggests 
that this effect is maximised at an optical depth of around 
k — 0.05, while at optical depths approaching zero the point 
masses should behave independently (in terms of microlens- 
ing), producing an effectiveness parameter of atot = 1- 

As the optical depth is increased, the population of 
caustics becomes more dense. This increase in density corre- 
sponds to an increase in the proportion of large light curve 
derivatives, as well as a decrease in the proportion of low 
derivatives. The increase in the contribution of fold caustics 
is reflected in both the increased effectiveness parameter, as 
well as the location of the derivative where the minimum KS 
difference is found. The relatively small size of this derivative 
highlights the fact that the total effectiveness parameter is 
a representation of microlensing flux variations of all levels. 

6.2 The effect of a mass function 

Previous work has demonstrated that the characteristic 
time-scale for variability scales as -J (m) , but is independent 
of the form of the microlens mass function (Witt, Kaiser 
& Refsdal 1993; Lewis & Irwin 1996). Also, the probabil- 
ity distribution of magnifications is a function only of the 
parameters of optical depth and shear (eg. Lewis & Irwin 
1995). This is demonstrated for the case of models 5 and 10 
in the left hand plot of Figure [|. We therefore expect that 
the distribution of microlensed light curve derivatives is a 
function of *J (m), but that atot is independent of the mass 
function. 



We investigate this by repeating some of the simulations 
of section ^ with a microlensing model of type SALPETER 
(models 9-13 in Table |l|). With the introduction of the mass 
function the mean mass of stars in the model has been low- 
ered. There is therefore an increase in the number density 
of stars in this model, and a corresponding increase in the 
number density of caustics. The increase in caustic density 
at constant optical depth is not expected to have an ef- 
fect on a to t because the rates of microlensing due to stellar 
proper motions and transverse velocity are increased equally. 
Figure ^| demonstrates this behaviour through comparison 
of derivative histograms that result from M SOLAR type 
simulations (solid lines) and SALPETER type simulations 
(dashed lines). The histograms of proper motion light curve 
derivatives are represented by the dark lines while the trans- 
verse motion histograms are represented by the light lines. 
The transverse velocities and the magnitudes of the proper 
motions that produce the microlensing in these histograms 
are each 300 km sec -1 . The histograms constructed from the 
SALPETER model show an increase in the number of large 
derivatives and a decrease in the number of small derivatives, 
demonstrating the increase in microlensing rate where the 
mean mass of stars is lower. The effect is applicable equally 
to both the proper motion histogram, and the transverse ve- 
locity histogram. This is demonstrated by the values of atot 
in TableH, which can be compared with those obtained in 
section [jJ~These values are also plotted in Figure |^ (thick 
light line). The results demonstrate that atot is independent 
of the mass function. 



6.3 The Effect of a Velocity dispersion 

Sections ^| and 6.2 have assumed the unrealistic, but easily 
interpreted case where all stars are given an equal speed in 
a random direction in the lens plane. These models produce 
distributions for light curve derivatives that result from the 
proper motion of stars, and from the equivalent transverse 
velocity that have KS differences of ~ 10 -2 . We need to 
check that this level of equivalence holds in the more realis- 
tic case where the stars are given a dispersion of velocities 
rather than a constant value. For models in this section, the 
components of the stellar velocities were each assigned at 
random according to a Gaussian distribution. We label this 
type of distribution by GAUSS. DISP. Simulations of type 
SALPETER were made at various optical depths which 
cover the range of interest in the case of Q2237+0305. 

With the inclusion of a velocity dispersion, a to t has been 
redefined as the equivalent transverse velocity divided by the 
average stellar speed. Table |^ shows the dependence of atot 
on k, as well as the minimised KS-difference and the deriva- 
tive at which it is found (in brackets). The trend, which is 
shown in Figure ^ (dashed line) is similar to the one obtained 
at the same optical depths in sections ^ and 3.2 . However the 
value of the total effectiveness parameter is larger at each op- 
tical depth considered. Kundic, Witt & Chang (1993) have 
shown that a velocity dispersion produces a higher average 
caustic velocity (with respect to the average stellar speed). 
The high velocity caustics produce a greater proportion of 
large light curve derivatives, and consequently leave a dearth 
of small derivatives. Importantly, the KS difference between 
the proper motion and scaled transverse histograms is again 
of order 10~ 2 in these models. The method outlined here 
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Figure 6. The probability density functions in both the cases of microlensing due to a transverse velocity (light lines), and due to 
stellar proper motions (dark lines). The probability density functions are shown for simulations of type M SOLAR (solid lines) and 
SALPETER (dashed lines). The optical depth in these simulations was k = 0.3, there was no applied shear (7 = 0). 



Table 4. Results from simulations of type M SOLAR, with ve- 
locity dispersions that are of type CONST. DISP. The optical 
depth in these simulations was k = 0.3. Details are found in the 
text. 



Shear 


Effectiveness 


KS 




Parameter (atot) 


Difference 


-0.4 


1.42±.05 


0.006 (.03) 


-0.2 


1.28±.04 


0.005 (.014) 


0.0 


1.08±.02 


0.012 (.007) 


0.20 


0.97±.02 


0.016 (.017) 


0.40 


0.92±.02 


0.009 (.029) 



can therefore be used to analyse the effect of proper mo- 
tions in realistic microlensing simulations of objects such as 
Q2237+0305. 



7 VARIATION WITH SHEAR 

The addition of an applied shear on the model starfield 
breaks the isotropy of the model. Kundic, Witt & Chang 
(1993) have shown that in the presence of an applied shear, 
the mean caustic velocity is dependent on the direction along 
which it is measured. In the case of microlensed flux varia- 
tion of a stationary source that results from stellar proper 
motions, the only component of the model that has a prefer- 
ential direction is the shear vector. The rate of microlensing 





x 
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-0.2 


APPLIED SHEAR 
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Figure 7. The relationship between the applied shear 7 and 
the effectiveness parameter atot- In these simulations, the optical 
depth was k = 0.3. 

produced through proper motions (by any measure) is there- 
fore not dependent on the direction in the source plane along 
which the light curve derivatives are sampled. The magnifi- 
cation distribution is also independent of the orientation of 
the shear (eg. Lewis & Irwin 1995), an independence that is 
demonstrated for models 14 and 17 in the right hand plot of 
Figure ^. However the rate of microlensing from the trans- 
verse motion of the source with respect to a static starfield 
that is perturbed by an applied shear is dependent on the 
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relative orientations of the transverse velocity and the shear 
(eg. Witt, Kaiser & Refsdal 1993; Lewis & Irwin 1996). The 
effect of an applied shear on the caustic network is to stretch 
the caustics along a direction that is at right angles to the 
shear (7 < 0). The result of this directional clustering is 
that a source moving in a direction that is parallel to the 
shear (7 > 0) experiences a higher microlensing rate. The 
directional dependence of the microlensing rate results in 
the effectiveness parameter also being dependent on the di- 
rection of the transverse motion. 

Table |i| shows the total effectiveness parameter at var- 
ious values of applied shear. The values were computed for 
a field having an optical depth of k = 0.3 in IMq stars 
(M SOLAR) . The stars in these models were given a proper 
motion of type CON ST.DI SP. The dependence of a to t on 
7 is plotted in Figure ^. This plot shows that the addition of 
a shear that is oriented at right angles to the source trajec- 
tory increases atot, while a parallel shear reduces a to t with 
respect to that of the same field in the absence of an applied 
shear. Table ^ also shows the minimised KS difference and 
the derivative at which it is found (in brackets). With the 
addition of an applied shear, the distributions of light curve 
derivatives resulting from proper motions, and from a static 
field in combination with the equivalent transverse velocity 
have best fit KS differences that are again of order 10~ 2 . 
Table ^ shows the variation of an me with shear. As in the 
case of variation with optical depth, an me follows a similar 
trend to that of atot- 

Figure M displays probability density functions of light 



Table 5. Results from simulations of type M SOLAR, with ve- 
locity dispersions that are of type CONST. DISP. The optical 
depth in these simulations was k = 0.3. Details are found in the 
text. 



Shear Effectiveness Parameter Effectiveness Parameter 
(« = 0.3) a HME (4#f > 0.15) a HME (iff > 1) 



-.40 1.36 1.30 

-0.20 1.27 1.24 

0.0 1.16 1.01 

0.20 1.02 0.98 

0.40 0.84 0.82 



curve derivatives produced both by a transverse velocity 
(light lines) and proper motions (dark lines) (each having 
a magnitude of 300 km sec) in the cases of 7 = —0.4 (solid 
lines) and 7 = +0.4 (dashed lines). As required the distribu- 
tion produced by proper motions is independent of the sign 
and therefore direction of the shear. However the transverse 
velocity produces a histogram of light curve derivatives that 
is dependent on the relative orientation of the source tra- 
jectory and shear. In the case where the shear is aligned 
at right angles to the source (7 < 0), lower derivatives are 
more common, and higher derivatives less common than in 
the case where the shear is aligned with the source trajectory 
(7 > 0). 
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Table 6. Values of the total effectiveness parameter for mi- 
crolensing models of type SALPETER corresponding to the im- 
ages of Q2237+0305. The models include a velocity dispersion of 
type GAUSS. DISP. 



Image 




atot 






degrees 


45 degrees 


90 degrees 


A,B (0.36,0.40) 


0.99±.08 


1.18±.10 


1.57±.10 


C (0.69,0.71) 


1.57±.05 


1.71±.07 


2.26±.27 


D (0.59,0.61) 


1.26±.05 


1.49±.12 


1.92±.15 



o 




20 40 60 80 100 

ANGLE (DEGREES) 

Figure 9. The relationship between the effectiveness parame- 
ter atot and the angle between the transverse velocity and shear 
for the microlensing parameters corresponding to the images in 
Q2237+0305. Images A/B, C and D are represented by the thin 
dark, dashed and thick light lines respectively. 

8 APPLICATION TO Q2237+0305 

As a final consideration we look at the values of atot corre- 
sponding to the parameters that describe the microlensing 
environment for the images of Q2237+0305. When mod- 
elling microlensing in Q2237+0305 we adopt the lensing 
galaxy model of Schmidt, Webster & Lewis (1998). This 
model includes the effect of the bar and produces the follow- 
ing microlensing parameters: image A (k = 0.36, |t*| = 0.40), 
image B (k = 0.36, |f| = 0.40), image C (« = 0.69, 
|7| = 0.71), image D (« = 0.59, \j\ = 0.61). Simulations 
of type SALPETER were made for each set of parameters 
(models 18-26 in Table Q) and a to t found in the cases where 
the angle between the transverse velocity and the shear is 0, 
45 and 90 degrees. The stellar proper motions were modelled 
by a distribution of type GAUSS. DISP. 

The results for atot are shown in Table |(| and plotted 
in Figure ^. As in section |^ these show that a to t is larger 
in the case where the transverse velocity is perpendicular to 
the source (7 < 0). The microlensing parameters (k and 7) 
increase in the order of images A/B, D and C. Figure [] there- 
fore outlines the general trend of atot which increases with 
optical depth and with a shear that is applied at right angles 
to the source trajectory. The results also demonstrate that 
for a given set of microlensing parameters, values of atot vary 
between the two extreme values for source trajectories that 
are oriented with respect to the shear at angles other than 
or 90 degrees. In the case of Q2237+0305, the four images 



are positioned approximately orthogonally with respect to 
the galactic centre. Therefore, the direction of galactic trans- 
verse motion with respect to the shear is approximately or- 
thogonal in images A and B to that in images C and D. An 
important point to be noted from Figure M is that when all 
images of Q2237+0305 are considered together, the relative 
rates of microlensing due to the transverse motion and stel- 
lar proper motions may be dependent on the direction of the 
galactic transverse motion. 



9 CONCLUSION 

We have combined the contouring algorithm of Lewis et al. 
(1993) with analytical expressions for both the derivative of 
the image magnification resulting from the proper motion 
of stars, and from the relative motion due to a transverse 
velocity of the galaxy. This combination has been used to 
construct probability density functions for point source light 
curve derivatives in both cases. We have compared the two 
classes of distributions by scaling the derivative in the trans- 
verse case to a value that minimises the KS difference. We 
find that the distribution of light curve derivatives where mi- 
crolensing is due to proper motion of stars is approximately 
equivalent to that of microlensing from a static lens with 
an appropriately scaled transverse velocity. The minimised 
KS differences are of order 10 -2 . The ratio of this appro- 
priately scaled transverse velocity to the mean speed of the 
stellar proper motions is defined as the total effectiveness 
parameter atot, and is a measure of the relative rates of mi- 
crolensing in the two cases. The value of atot is dependent 
on the microlensing parameters (« and 7), but is indepen- 
dent of the assumed microlens mass function. We find that 
a to t increases with optical depth in the absence of a shear, 
increases with shear that is aligned at right angles to the 
transverse source trajectory, and decreases with a shear that 
is aligned along the transverse source trajectory. The latter 
results are particularly important in the case of Q2237+0305 
where 4 images are placed orthogonally with respect to the 
centre of the lensing galaxy and so have varying directions 
of galactic shear with respect to the transverse motion. It 
may therefore be possible to constrain the direction of trans- 
verse galactic motion (assuming k and 7) from the observed 
microlensing rate. 
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APPENDIX A: EXPRESSIONS FOR THE DERIVATIVES 

We have defined the following quantities for use in the evaluation of ^jf-) pro per and ^jf-)tran 

j (x 2 - X 3 2 ) 2 - (Xl - x{) 2 



A 3 = 
B J = 2m 



j (a?a - x 3 2 )(xi - x\) 

((x 2 ~ x 2 y + ( Xl - x{yy 

C ] = 2m i (x 1 -xi) 



D J 



A 



c 



2m 3 {x 2 — x 2 ) 



3(x 2 - x\f - (xi - x{) 2 

■ {(x 2 - x ] 2 y + (xi - x{yy 

3(x 1 - x{) 2 ^ (x 2 - xj) 2 

(( Xl - x{) 2 + (x 2 ~ xiyy 



J2 AJ B = Y1 B3 

3=0 3=0 
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£V D^D* 



3=0 



3=0 



In terms of these quantities, the full expression for the derivative -jf-) P roper is: 



~dt 



C( 7 + A) + BD 



-' 2 "r [ B2 1' {1 + m x E (ft 1 + (T + A )) AJ + BB3 ) < + ( BAJ " C 1 + (7+ A)) B 3 ) n 
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/ N * \ 
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AT, 

+2 M 2 J2 (° J (7 + A) + BD 3 ) v{ 
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The definitions Al can also be used in the expression for the static field light curve derivative: 

mB-(l + ( 7 + /l)) > 



^f) = 24((D( 1 + A)-BC)x 

Uib J tran 



B - (1 - ( 7 + A)m) 
B 2 - (1 - (7 + AY) 



+ (C( 7 + A) - BD) x 



B 2 -(l-(7 + A) 2 ) 



x v u 
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